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We study a disordered classical Heisenberg magnet with uniformly antiferromagnetic interactions 
which are frustrated on account of their long-range Coulomb form, i.e. J{r) ~ — Alnr in d = 2 
and J(r) ~ A/r in d = 3. This arises naturally as the T —>■ 0 limit of the emergent interactions 
between vacancy-induced degrees of freedom in a class of diluted Coulomb spin liquids (including the 
classical Heisenberg antiferromagnets on checkerboard, SCGO and pyrochlore lattices) and presents 
a novel variant of a disordered long-range spin Hamiltonian. Using detailed analytical and numerical 
studies we establish that this model exhibits a very broad paramagnetic regime that extends to very 
large values of A in both d = 2 and d — 3. In d = 2, using the lattice-Green function based 
finite-size regularization of the Coulomb potential (which corresponds naturally to the underlying 
low-temperature limit of the emergent interactions between orphan-spins), we only find evidence 
that freezing into a glassy state occurs in the limit of strong coupling, A = oo, while no such 
transition seems to exist at all in d = 3. We also demonstrate the presence and importance of 
screening for such a magnet. We analyse the spectrum of the Euclidean random matrices describing 
a Gaussian version of this problem, and identify a corresponding quantum mechanical scattering 
problem. 

PACS numbers: xx 


I. INTRODUCTION 

The appearance of novel magnetic phasesip— generally 
contains as one ingredient the ability of the system to 
avoid conventional (semi-)classical ordering. In this con¬ 
nection, the role of several factors has been extensively 
explored. These include low dimensionality and the re¬ 
sulting enhancement in the effects of quantum and en- 
tropic fluctuations, geometrical frustration, whereby the 
leading antiferromagnetic interactions compete with each 
other on lattices such as the kagome and pyrochlore lat¬ 
tice, and the presence of quenched disorder, which dis¬ 
rupts any residual tendency to conventional long-range 
order. Each of these has given rise to research efforts 
spanning decades of work. 

Here, we study a model with a new combination of 
some of these ingredients. The focus of our study is 
a disordered classical Heisenberg magnet with antifer¬ 
romagnetic interactions which are frustrated on account 
of their long-range Coulomb form at long-distances, i.e. 
J{r) ^ —Aln(r/£) in d = 2 (where £ is a length-scale 
of order the system-size) and J{r) ^ A/r in d = 3. This 
Coulomb form of the Heisenberg couplings arises natu¬ 
rally as the r —>■ 0 limit of the emergent entropic ex¬ 
change interactions^ between vacancy-induced “orphan- 
spin” degrees of freedom^^— in diluted Coulomb spin liq¬ 
uids, and presents a novel variant of a disordered long- 
range spin Hamiltonian with connections to Euclidean 
random matrices. The coupling constant A is determined 
in any given system by the microscopic details of the un¬ 
derlying Coulomb-spin liquid, while the spin degrees of 


freedom in the model we study are related to the physi¬ 
cal orphan-spins of the underlying diluted magnet. Our 
focus here is on studying the range of behaviours possible 
in the T —>■ 0 limit by mapping out the phase diagram of 
our Coulomb antiferomagnet as a function of A. Erustra- 
tion arises naturally in the model under consideration, as 
any triplet of spins are mutually coupled antiferromag- 
netically but without the randomness in sign of, say, the 
Sherrington-Kirkpatrick model^. Also, unlike the latter 
case, the interactions are long-ranged but not indepen¬ 
dent of distance. 

Our motivations for studying it include having been 
led to this model in a previous investigationiS of diluted 
frustrated magnets exhibiting a Coulomb spin liquid at 
low temperature. The model is in this sense natural, 
appearing as the zero-temperature limit of a disordered 
frustrated magnet. The corresponding experiments are 
on the material known as SCGO, which triggered the in¬ 
terest in what we now call highly frustated magnetism in 
the late Its behaviour at very low temperatures is 

still not very well understood, e.g. the observed glassiness 
even at very low impurity densitie o^^d^ , which appears 
to involve only the freezing of a fraction of its degrees of 
freedom. We will return to this point in Sec. lVH 51 While 
exhibiting a classical Coulomb spin liquid regime, the dis¬ 
order in this system leads to the emergence of new, frac¬ 
tionalized, degrees of freedom, the so-called Orphans^i^, 
which interact via an effective entropic long range inter¬ 
action mediated by their host spin liquid^. 

We believe that as such, it can be of interest as a 
generic instance of the interplay of strong interactions 




2 


and disorder in magnetism. In particular, it develops the 
strand of thought of how disorder in a topological sys¬ 
tem characterised by an emergent gauge field can nucle¬ 
ate gauge-charged defects, with the pristine bulk medi¬ 
ating an effective interaction between them. Long-range 
Coulomb interactions like the one studied here are then 
as natural as the algebraically decaying RKKY interac¬ 
tions in metallic spin glasses. 

Our central results are the following. First we use the 
results of previous worb^, to work out in detail the key 
features of this T —>■ 0 limit, and demonstrate that this 
limit is characterized by a single coupling constant A, 
which is, in principle, determined by the geometry of 
the underlying spin-liquid. Second, our extensive Monte 
Carlo simulations for d = 2 reveal no sign of any freezing 
or ordering transition up to very large coupling strengths. 
At the same time, within a self-consistent Gaussian ap¬ 
proximation, we find that there does appear such a tran¬ 
sition at infinite coupling in d = 2 but not in d = 3. 
This transition is very tenuous, in that it is replaced by 
a more conventional ordering transition in a finite system 
depending on the choice of how to regularize this long- 
range interaction in a finite lattice: the finite-size lattice 
regularization that is most natural from the point of view 
of the T —>■ 0 limit of the underlying diluted magnet gives 
rise to freezing into a glassy state at A~^ = 0, while other 
regularizations replace this glassy state by a conventional 
ordering pattern. The Coulomb antiferromagnet there¬ 
fore remains highly susceptible to perturbations, just like 
many other frustrated magnets^. 

We also study the spectrum of the interaction matrix 
of this random Coulomb antiferromagnet, which provides 
an instance of an Euclidean random matri x^^d^ , in that 
its entries are obtained as a distance function between 
randomly chosen location vectors. We find two qualita¬ 
tively distinct regimes. On one hand, at low energies in 
the low-density limit, eigenfunctions are localised, with 
the lowest energy states as pairs of neighbouring spins 
the probability distribution of which we compute. Be¬ 
yond this extreme low-density limit, more complex lat¬ 
tice animals appear in this regime. On the other hand, at 
high energies, the modes correspond to long-wavelength 
charge density variations with superextensive energy. In 
between, we find no clear signature of a well-defined mo¬ 
bility edge in this Coulomb system. 

Another interesting aspect of the uniformly antiferro¬ 
magnetic interactions is that they permit a variant of 
screening to appear in this Coulomb magnet, which has 
no correspondence with other long-range magnets such as 
the Sherrington-Kirkpatrick model. Our analysis of this 
screening further leads us to an identification of the cor¬ 
relations of the random Coulomb antiferromagnet with 
the properties of the zero-energy eigenstate of a quan¬ 
tum particle in a box with randomly placed scatterers. 

Returning to experiments, we note that the uniform 
magnetic susceptibility of SCGO will of course be dom¬ 
inated by the Gurie tail 1/T) produced by these or¬ 
phan spins at low temperature. Both in d = 2 and d = 3, 


the full susceptibility, when vacancies are placed at ran¬ 
dom, is that of independent orphans to a good approx¬ 
imation despite the long-ranged interaction present be¬ 
tween them. This persists down to the lowest temper¬ 
atures not only because of the screening of the interac¬ 
tions at hnite physical temperature, and because the size 
of the Goulomb coupling derived from the entropic inter¬ 
action is comparatively weak, but also because the phys¬ 
ical orphan spins are related to the degrees of freedom in 
the Goulomb antiferromagnet via a sublattice-dependent 
staggering transformation, so that the uniform suscep¬ 
tibility of the physical orphan-spins corresponds to the 
staggered susceptibility of the degrees of freedom of our 
Goulomb antiferromagnet, and therefore remains largely 
unaffected by the fact that the total (vector) gauge charge 
of our Goulomb antiferromagnet vanishes. 

The remainder of this paper is structured as follows. In 
Section |TT1 we first provide a self-contained review of ear¬ 
lier work on vacancy-induced effective spins in a class of 
classical antiferromagnets on lattices consisting of “cor¬ 
ner sharing units”, and then build on this to provide a 
careful derivation of the T —>■ 0 limit of the emergent en¬ 
tropic interactions between orphan spins and use this to 
dehne our model Goulomb antiferromagnet. After outlin¬ 
ing our analytical and numerical approaches in Sec. IIIII 
we present the results obtained in d = 2 and d = 3 
(Sec. IIVI) . Sec. IVl contains the analysis of the problem 
in terms of a Euclidean random matrix while the role of 
screening and the connection to a scattering problem are 
discussed in Sec. ED We conclude with a discussion of 
these results, and relegate sundry details (such as dicus- 
sions of the fully occupied lattice and the ordered state 
seeded by a certain finite-lattice regularization of the two- 
dimensional Goulomb interaction) to Appendices. 


II. THE RANDOM COULOMB 
ANTIFERROMAGNETIC HAMILTONIAN 

We thus study a classical Heisenberg model 

id = — JijUi ■ Uj. (1) 

ij 

where Jij takes on a Coulomb form, 

Jij = -A\og{rij/C) (d=2) (2) 

= A/r,, (d = 3). (3) 

This form with C larger than any Vij has the property 
that the interactions are uniformly antiferromagnetic as 
well as long-ranged. 

We need to supplement this by defining the degrees of 
freedom, unit vectors rii, appearing in Eq. [TJ We con¬ 
centrate on the case where their locations, denoted by i 
are chosen randomly on a square (cubic) lattice in d = 2 
(d = 3), at a dimensionless density of x spins per lattice 
site. 
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FIG. 1: Illustration of the Orphan spin arising from the in¬ 
troduction of non-magnetic impurities. Its effective moment 
is half of that of a free spin. 


For long-range interactions like this Coulomb inter¬ 
action, choices about boundary conditions or ensemble 
constraints can be considerably less innocuous than for 
short-range systems. In order to illustrate this, and to 
make natural choices for these items, as well as for moti¬ 
vation of our study, we discuss the derivation of a random 
Coulomb antiferromagnet as an effective Hamiltonian of 
a diluted Coulomb spin liquid next. 


These local constraints lead to an effective description in 
terms of a theory of emergent electric fields that obey a 
Gauss law. To see this, we define electric fields E2 = tj-SS 

on links I, where is a spatial unit vector that points 
from the A- to the H-sublattice of the premedial lattice 
end of this link. The ground-state condition then trans¬ 
lates to the statement that the lattice-divergence of this 
electric field vanishes at each site K for each a. The key 
idea of this effective description is that the coarse-grained 
(entropic) free energy density depends quadratically on 
the local electric field, and deviations from the vanishing 
divergence condition amount to the appearance of vector 
Coulomb charges^. These emergent gauge charges are de¬ 
fined for each lattice point Kl of the bipartite premedial 
lattice: 


Qki = vi^) X! (6) 

feKl 

and the staggering factor, r;(Kl) = -1-1 if is an A- 
sublattice site of the premedial graph and —1 otherwise. 
Since each microscopic spin contributes with opposite 
signs to the vector charge on two neighbouring simplices, 
the total gauge charge of a system without boundaries 
must vanish in every configuration of the system 

= ( 7 ) 


A. Orphan spins and their interactions in diluted 
Heisenberg antiferromagnets 

We thus begin by providing a self-contained review of 
earlier work on vacancy-induced effective spins in a class 
of classical frustrated antiferromagnets on lattices con¬ 
sisting of “corner sharing units”. The centers of these in 
turn define a so-called premedial lattice, which is bipar¬ 
tite in practically all instances of popularly studied clas¬ 
sical Heisenberg spin liquids^. A simple model of nearest 
neighbour antiferromagnetically interacting spins on such 
lattices can be written as 

^ feKI 

where the summation in the alternate form of the Hamil¬ 
tonian is carried over the corner sharing simplices K, 
which might be tetrahedra, as e.g., in a pyrochlore lattice, 
triangles in a Kagome lattice, or a combination of both 
as in the case of SCGO, and the spins of the frustrated 
magnet are now labeled by /, the links of the bipartite 
pre-medial graph (whose sites correspond to the centers 
of the simplices of the original lattice, and links I cor¬ 
respond to sites of the original lattice). When written in 
this form, it is clear that ground-states are characterized 
by the constraints: 

E‘^r=o,v^, (5) 


Kl 

This very natural condition-akin to the charge-neutrality 
of the full universe, and in our case unavoidable due to 
the microscopic origin of the emergent gauge charge - will 
be explicitly imposed in our Monte Carlo simulations of 
the system. 

The mapping of the pure system to an emerging gauge 
field theory at low temperatures makes clear that gener¬ 
alized “vector charges”, Qk , are generated thermally as 
a consequence of the violation of the ground state con¬ 
straints. The constraint Eq. [S] is also unavoidably vio¬ 
lated in the presence of non-magnetic impurities (Fig. [T]) 
whenever all but one spin of a given simplex are sub¬ 
stituted for by vacancies (simplices containing at least 
two spins can in general satisfy the zero total spin condi¬ 
tion and such simplices do not host a vector charge in the 
T —>■ 0 limit). Indeed, when all spins but one in a simplex 
are replaced by vacancies, the result is a paramagnetic 
Curie-like response^i^d^, which dominates the suscepti¬ 
bility response at low temperatures. The lone spins on 
these defective simplices, which serve as the epicenter of 
this paramagnetic response, were baptized Orphans (Ref. 
d) in the first studies of this effect. 

The field theory developed in Refs. M extends the 
self consistent gaussian approximation (SCGA)i^, a the¬ 
ory successful in describing low temperature correlations 
on the undiluted systems, to incorporate the effects of 
dilution and study the physics of these orphan spins at 
non-zero temperature in a manner that treats entropic 
effects on an equal footing with energetic considerations. 
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In its original form the SCGA replaces the hard con¬ 
straint on the spins norm, Sf = S'^, by the relaxed soft 
spin condition on their thermal average (Sf) = S^. The 
key insight of Refs. that led to the detailed analyti¬ 
cal understanding summarized below, was the following: 
While it is sufficient to treat in this self-consistent Gaus¬ 
sian manner all spins other than the lone orphan spin in 
a simplex in which all but one spin has been replaced 
by vacancies, this is much too crude an approximation 
for the orphan-spin itself, which must be treated with¬ 
out approximation as a hard-spin obeying 
Remarkably, the resulting hybrid field theory continues 
to be analytically tractable when the number of orphan 
spins is smaMd^. With just one orphan present in a sam¬ 
ple with an external magnetic field of strength B along 
the 0 axis, the theory predicts that this orphan spin sees 
a magnetic field 5/2, with the other half of the exter¬ 
nal field screened out by the coupling to the bulk spin- 
liquid. The resulting polarization of the orphan serves as 
a source for an oscillating texture that spreads through 
the bulk. The net spin carried by the texture cancels 
half the spin polarization of the orphan-spin, resulting 
in an impurity susceptibility corresponding to a classi¬ 
cal spin S/2. With more than one orphan present, the 
spin-textures seeded by each orphan mediate an effective 
entropic interaction between each pair of orphan spins. 

The effective action for a pair of orphans is predicted 
in this manner to have the form 


-/3Jeff(r,r)ni • n 2 , (8) 


where ft are unit-vectors corresponding to the directions 
of the orphan-spins in a given configuration. The ex¬ 
change coupling has a particularly simple form in the 
large separation limit 


PJeS 


-r]{ri)ri{r2) 


{Qm{ri) ■ Q^{r2)) 
(Qki ■ 


(9) 


which involves only “charge-charge” correlations calcu¬ 
lated in the pure system'. 




• {fi - rs)) 




( 10 ) 


The denominator of Eq. ®, behaves at low temperatures 
as • <5^) = T/ J from equipartition. 

For orphans in d = 2, one finds: 


the entropic screening length ^ent = 1 /k ^ ^/'JT sepa¬ 
rates two regimes, algebraic /C(r) ^ —1/r and screened 
/C(r) ~ exp(—Kr). 

In the physical system, at any nonzero temperature, 
this is thus a ‘short-ranged’ interaction on account of the 
finite screening length which, however, diverges as 1/\/T. 
In this article, we are interested in the limit of T = 0, 
where the interaction takes on the novel - for magnetic 
systems - long-range Goulomb form. 


B. Model Hamiltonian 

In the limit of T —>■ 0, we are thus led by these con¬ 
siderations to Goulomb interactions between the vector 
orphan spins, which we here study in detail. For simplic¬ 
ity, we consider unit-vector spins n at random locations 
in a periodic hypercubic lattice of linear size L with oc¬ 
cupancy probability x, corresponding to an underlying 
spin liquid on the checkerboard and “octochlore” lattices 
of corner-sharing units involving 2'^ spins in d dimension. 

In what follows, we will get rid of the sublattice factors 
that affect the sign of the effective interaction by invert¬ 
ing all unit-vectors placed on the B sublattice. In other 
words, we identify Sfii with ryi^orphan,*, where S'orphan.i 
is the orphan spin on the simplex labeled by i in the 
underlying diluted frustrated magnet. 

This gives us a “random Goulomb antiferromagnet” in 
which unit-vector spins interact with an exchange cou¬ 
pling that is always antiferromagnetic but of a long-range 
Goulomb form at large distances. For a classical system, 
this transformation is innocuous, but note that natural 
observables cease to be so under this mapping - e.g. the 
orphan spin contribution to the uniform susceptibility of 
the underlying diluted magnet is now given by the stag¬ 
gered susceptibility of our Goulomb antiferromagnet. 

As is usual for entropic interactions in the limit of 
T —>■ 0, the strength of their coupling. A, is fixed by the 
microscopies of the model from which they have emerged. 
In this work, we are interested in exploring the generic 
behaviour of such models - in particular, identify possi¬ 
ble phases - and thus allow the coupling A to be vari¬ 
able. For completeness, we mention that A = ^ for the 
checkerboard lattice. 

This therefore leads to the form of H at the begin¬ 
ning of this section, Eq. [TJ To make Eq. [2] dimensionally 
unambiguous we write: 


Jy = -Alog(r-y/£) (d = 2) 


Jes{ri -r2,T) = r]{ri)ri{r2)TJ{\ri -r2\/fent) (H) 

with an entropic screening length ^ent = 1 /k ~ I/VT 
separating two regimes for J{Kr). For /tr <C 1 a logarith¬ 
mic one, ~ — log(K7’); and for Kr 3> 1 a screened 

regime, J{nr) ~ exp(—Kr). Analogously in d = 3, 

JeSin - f2,T) = ?7(ri)?7(r2)r^/^/C(|ri - f2|/^ent) (12) 


with C conveniently set to a value of order the system size 
L so that Jij > 0 always. In the above language, with 
sublattice factors ij absorbed into the definitions of Ui, 
the zero gauge-charge constraint imposed by the micro¬ 
scopic origin of this effective model now translates to the 
constraint that Uj = 0 in every allowed configuration 
of our Goulomb antiferromagnet. This constraint in fact 
can also be imposed by adding an infinitely strong inter- 
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action acting equally between all spins. This equivalence 
renders the detailed choice of C immaterial. 

We note an interesting scale-invariance of this model in 
the limit of small densities of spins. This scale invariance 
is inherited from that of the logarithmic function under 
scaling transformations: J{nr) = log(K) + J{r), together 
with net charge neutrality Eq. ©: 

^n, = 0, (13) 

i 

implies that the extra term log(K) gives a temperature- 
independent contribution to the action determined by 
= —N/2. The partition function thus 
only picks up a constant factor: 

log{K.) N / 2 ^ ^ 

It also means that, rather unusually, in the continuum 
limit X —>■ 0 the partition function is a scaling function 
depending on the randomly chosen orphan locations only 
scaled by their mean separation. Lattice discretisation 
effects at finite x break this equivalence. The scaling 
transformation for the model in three dimensions gives 
J{Kr) = J(r)/K^ what implies for the partition function 
a rescaling of /?: 

Z'{P) = ^ Z{P/k). (15) 

For Coulomb interactions in a finite-size system, var¬ 
ious choices of the interaction yield the same large- 
distance form in the limit L —> oo. The most natural 
form from the point of view of the effective field the¬ 
ory predictions for emergent interactions between orphan 
spins is the Fourier transform of the inverse of the lattice 

d 

Laplacian, d—^ coski'. 

i=l 


J{xij) 




d 

d — cos ki 

i=l 


(16) 


This we call the lattice Green function (LGF), and our 
most detailed studies are carried out with this form of 
the interaction. 

Alternatively, one can work directly with the Goulomb 
form, e.g. for d = 2: 


J(ry) =-log . (17) 

with C = L/-\/2. This form agrees with the LGF inter¬ 
actions at short distances (see Fig. El). 

The issue of how to impose the boundary conditions, 
and therefore how to compute r^ , turns out to make 
much difference on the results for a finite system, as we 
shall see. The choices of either 

Tij = \n - gI = yjxl+yjj, (18) 



FIG. 2: J(x,y) used in the simulations in d = 2. 


with Xij = min(|xi — Xj\,L — \xi — Xj\) or 

n, = + sin» 

(19) 

result in different behavior for the system, which will be 
explained in more detail in the results section. We refer 
to these choices as periodised, and smoothed, logarithms, 
respectively. The latter is very close to the LGF, while 
the former maintains a finite difference to it at the pe¬ 
riodic boundary, where it is not differentiable for any L 
(Fig. [2|). It is easily seen why this finite difference is in¬ 
dependent of L, if one compares the smoothed log to the 
periodized Log, approximatelly equivalent to comparing 
the LGF with the Log. Looking, e.g., at the midpoint of 
one edge (x^ = L/2, yij = 0) one finds: 

( jLGF _ jLog)(L/2, 0) log(7r/2), (20) 

where the subindex L emphasizes that we are looking at 
the respective forms of the interactions in a finite system 
of size L. 

Note, again, that adding a constant to the interaction 
(in d = 2), e.g., by changing the denominator of Eq. (fTTl) . 
leaves the interaction unchanged due to the global charge 
neutrality constraint. 


III. METHODS 

The analysis of spin systems with the potential for 
glassy phases is a delicate endeavour as equilibration 
of large systems is elusive. Existence and determina¬ 
tion of a transition temperature is usually a controversial 
issueiid^. Since our system has long ranged interactions, 
boundary effects can cause yet more trouble. This is why 
we combine analytical with numerical methods, as well 
as mappings to other problems which have received at¬ 
tention in a different context previously. 

Numerically, we study the behaviour of this model 
through Monte Garlo (MC) simulations, and analyti¬ 
cally in the self-consistent Gaussian (“large-m”, also de- 
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noted in the following as LM approachi^Ti^i) approxima¬ 
tion, where the parameter A mimics an inverse tempera¬ 
ture. Our MC simulations directly impose the constraint, 
Eq. ([7]). For that we initialize the system in a random 
configuration of vanishing total spin, and the update 
movements on the system consist of selecting an arbi¬ 
trary pair of spins, and rotating them around the axis 
determined by their vectorial sum. A MC simulation of 
the same system with strictly positive interactions, with¬ 
out this constraint on the total spin has been also inves¬ 
tigated, and the conclusion is that while the relaxation 
time increases, the system still prefers to stay close to 
the manifold of vanishing total spin. 

The LM approach consists of considering spins with m 
components and letting to —>■ oo. This is formally equiv¬ 
alent to the soft spin approximation and it only gives in 
principle information about the infinite number of com¬ 
ponents limit, but this can be understood as the 1st term 
in an expansion of the 0(m) model. It has been very suc¬ 
cessful in the analytical study of correlations in highly 
frustrated spin systems!^, being able to reproduce the 
main features of the on-going phenomena, such as exis¬ 
tence of long range dipolar correlations at T = 0, charac¬ 
terized by the presence of “pinch points” in the structure 
factor—. 

The LM approach allows an analysis of the system both 
at finite coupling strengths A < oo, and at A = oo. The 
study of glassiness with this approach has been already 
undertaken in a variety of model o^^’^^ , and we will be 
following a similar methodology. Correlations are com¬ 
puted through the matrix: 


A fundamental quantity at A = oo within the LM 
approach is the number of zero eigenvalues, toq, of the 
matrix Bij = Jij -\- hidij; it can be showiti^ that the 
ground state spin vectors span an toq dimensional space. 
This quantity should scale with the number of particles 
in the system as toq Furthermore, as was shown 

in Ref. the same exponent controls the scaling of the 
spin glass susceptibility for the ground state configura¬ 
tion: XSG ^ N^~^. 

The main quantity of interest in our study will be 
the spin glass susceptibility (square brackets here and 
throughout indicate the disorder average). 


Xscik) 


1 

N 


• fij)'^ cos k ■ (r) 



(24) 


obtained in the MC simulations through the overlap ten¬ 
sor—: 



1 

N 




(25) 


where greek indices refer to the spin components, while 
the indices 1, 2 refer to two independent replicas of a dis¬ 
order realisation. This might be interpreted as the over¬ 
lap of a spin configuration with itself after an infinitely 
long time. Since the onset of glassiness can be also un¬ 
derstood as a divergence of the equilibration time, the 
nonvanishing of this order parameter signalizes the tran¬ 
sition. 

The spin glass susceptibility in terms of this tensor is: 


B^j — iJ-ij hi6ijj (^1) 

and are given by: 

Cij = (22) 

These can be computed once the Lagrange multipliers, 
hi, are determined through the set of nonlinear equations: 


Xscik) = 




Q 


q:,/3 


(26) 


We follow the usual practice to determine the spin glass 
transition by computing a finite system correlation length 
associated to the susceptibility above. The Ornstein- 
Zernike form for correlations gives: 


C,, = 1. (23) 

For comparison between LM and MC, we scale observ¬ 
ables and couplings with to so that their small-coupling 
(“high-temperature”) forms agree. 

The point A = oo is treated within the LM approach 
by determining the (uniquei^) ground state through a 
local field quench algorithm^!. This algorithm is based 
on the fact that if the number of spin components, to, 
is large enough (larger than i/ 2N —) , then a system of 
spins with to components is effectively equivalent to the 
corresponding system in the limit to —>■ oo. The algo¬ 
rithm then consists of taking a system of N spins with 
TO > \/2N components initially randomly oriented, and 
then iteratively aligning each spin with its local field. 
This procedure is expected to converge to the unique 
ground state, from which all the quantities of interest 
can be computed. 


( \ 1/2 

Xsg(O) — 1 I (27) 

XScikmin) / 

and near the transition, the finite size scaling prediction 
is expected to be: 

|i = A(Li/'^(l/A-l/A,)), (28) 

while the susceptibility should follow: 

|^=r(Li/‘^(l/A-l/A,)), (29) 

Notice that these scaling relations only hold if there 
exists a crossing of finite size correlation length curves 
for different system sizes at an unique finite coupling 
strength value. The absence of such a crossing at a finite 
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FIG. 3: Spin glass susceptibility (top) and correlation length 
(bottom) for the LGF interaction. Several system sizes are 
indicated by different colours. Circles indicate (error bars 
are of the order of the circles size) MC simulations, while 
lines are from the LM approach-correlations are stronger for 
Heisenberg spins than the ’soft’ LM spins throughout. The 
insets show scaling collapse for LM for IjA^ = 0. 

Ac indicates the absence of a phase transition. Nonethe¬ 
less a phase transition at Ac = oo cannot thus be ruled 
out and the LM approach allows an analysis in this situ¬ 
ation. The scaling relations predicted to hold in this case 
(Ac —>■ oo) are: 

= lA), = (30) 

The exponent ^ here is the one previously introduced 
for the scaling of the number of zero eigenvalues of the 
matrix B with the number of particles in the system. 


IV. RESULTS 
A. Two dimensions 

The two approaches (MC and LM) yield a broadly con¬ 
sistent picture for each of the interactions studied. We 
conduct an analysis of a possible freezing transition in the 


FIG. 4: Scaling of the number of zero eigenvalues (mo) of the 
matrix B defined in the text and of the spin glass susceptibil¬ 
ity (insets) with the number of particles for the LGF (top), 
and Log (bottom) interactions. 


model by measuring the spin glass susceptibility and try¬ 
ing to identify the transition through a finite size scaling 
of its associated correlation length. Other observables 
such as the specific heat or the uniform susceptibility 
were also studied, though these do not indicate any of 
the conventional orderings. 

The results from MC simulations and LM calculations 
are shown on Fig. [3] for the system with LGF as interac¬ 
tion for a fixed density x = 0.10 of particles. In each case 
the number of disorder realisations simulated was 200. 

Globally, correlations are stronger for the MC simula¬ 
tions on Heisenberg spins compared to the LM results. 
This is in keeping with the general lore that a lower num¬ 
ber of spin components is conducive to spin freezing, as is 
well known from the comparison of Ising and Heisenberg 
spins. 

In the broad range of coupling strengths considered by 
our analysis, no unique crossing for the different system 
sizes of the correlation length curves can be identified. 

The LM analysis at A = oo yields the exponent /i as 
indicated in Fig. |4l This seems to have the same value, 
/j, Ri 0.3 for both the LGF and Log interactions. 

The exponent value fj, = 0.3 is used as input, together 





































FIG. 5: Spin glass susceptibility (top) and correlation length 
(bottom) for the Log interaction, as computed on the MC 
simulations (points) or with the LM approach (lines). The 
insets show the corresponding scaling collapses. 


with the assumption that Ac = oo for the LGF, in at¬ 
tempting a scaling collapse of the LM data. The expo¬ 
nent ly was determined by a fitting procedure with the 
scaling relation, Eq. dsni), only using data for the corre¬ 
lation length. The resulting scaling collapse is shown on 
the inset of the lower panel of Fig. [3l where v = 0.68(1) 
is obtained. Finally, we use all these exponents on the 
predicted scaling relation for the susceptibility (the re¬ 
sult is shown on the inset of the upper panel of Fig. [3|). 
The available data from the LM calculations indicates 
therefore a freezing transition at = oo for the diluted 
model with LGF as interaction in two dimensions. 

The Log interaction turns out leads to a dramati¬ 
cally differing behaviour! This is a surprising result, as 
the interactions only differ appreciably at large distances 
(Fig. [2|). Fig. [5] shows the results for the observables of 
interest as obtained from MG simulations and LM cal¬ 
culations, respectively. Here again we fix the density of 
particles x = 0.1, and consider 200 disorder realisations. 
A clear crossing of the correlation length curves for differ¬ 
ent system sizes occurs and scaling collapses of the data 
are possible, which are shown together with the corre¬ 
sponding critical exponents as insets. 
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FIG. 6: Disorder-averaged pair correlations with a spin at the 
origin as a function of relative coordinates, centers of circles 
indicate the position of the spin, its radius gives the mag¬ 
nitude, with red (black) denoting positive (negative) correla¬ 
tions. The central red circle thus reflects {Si) — 1. The upper 
panel indicates the result for the LGF with A = 100, while 
the lower corresponds to the Log with A = 20. Data shown 
from MC is in agreement with LM (not shown). 


To study more closely this effect, we consider the pair 
correlations as a function of the relative coordinates of 
the pairs, averaged over disorder realisations (Fig. [51). 
The profile is isotropic for the LGF with only the 1st few 
nearest neighbors significantly antiferromagnetically cor¬ 
related. On the other hand, the Log interaction yields 
strongly anisotropic behavior (the interaction itself is 
anisotropic) and this seems to be responsible for what 
we see as a “glassy phase transition” emerging from the 
“splaying out” of the susceptibility curves. The absence 
of glassiness is explained in more detail on Appendix 
where we expose how the pair correlation profile helps us 
in defining an appropriate susceptibility for the case at 
hand, which is shown to diverge in the thermodynamic 
limit. It turns out that this reflects not the existence of 
true glassiness but a transition closer to conventional or¬ 
dering. Note that the gross features of the correlations 
(Fig. |5| lower panel) follow if one frustrates the pairs at 
the kink (Fig. |2|) of the Log interaction, which form a 
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frame at half the system size. The set of points which in 
turn are on the “frames” of 0{L) points on the first frame 
yield the cross shaped set of ferromagnetically correlated 
sites centred on the origin. 

Note that such finite-size differences appear to be ab¬ 
sent in previous studies in d = 1—; they appear to 
be a consequence of the anisotropic nature of our peri- 
odised Log interaction with its non-analytic minimum at 
maximum separation. By contrast, the “smoothed Log” 
(Fig. [2]) that also respects the periodic boundary condi¬ 
tions essentially reproduces the LGF interaction results. 


1. The fully covered square lattice 

For completeness, we have also analysed the situation 
for a fully occupied lattice. In this case we observe that 
the LGF interaction leads to conventional (Neel) anti¬ 
ferromagnetic order, while the Log leads to a “striped” 
phase. This can be understood from a theorem in Ref. [2^ 
which states that the ground state of the system is de¬ 
termined by the minimum of the Fourier transform of 
the interaction. This is explained in more detail on Ap¬ 
pendix |B] 


B. Three dimensions 

We analyse the diluted cubic lattice considering a den¬ 
sity of particles x = 0.0625, and again considering the 
model Hamiltonian of Eq. m, with interactions now re¬ 
stricted to be the LGF as given by Eq. (1161) . Both Monte 
Garlo simulations and LM calculations cover several sys¬ 
tem sizes with 100 distinct disorder realisations each. 
The main focus is on the possibility of a glassy phase 
and the spin glass susceptibility and corresponding cor¬ 
relation length are computed. Our prior discussion of the 
finite size scaling relations still holds, and one determines 
the transition as an unique crossing of the finite size cor¬ 
relation length curves. Instead of this we observe (Eig.[7]) 
a trend for the crossings to shift towards larger values of 
A as the system size increases, similar to the situation in 
two dimensions. 

No good scaling collapse was obtained. A freezing tran¬ 
sition in this system at a finite coupling strength therefore 
appears unlikely, though a more careful finite size scaling 
analysis of the crossings is necessary to give a definitive 
answer. 

A LM study at A = oo reveals that the exponent for 
the scaling of zero eigenvalues of the matrix B with sys¬ 
tem size yields fi = 0.33, in agreement with the predic¬ 
tion in 3 dimensions for a short ranged interacting sys- 
tem^i. Using of this exponent and the scaling relations 
at A = oo does not lead to a good scaling collapse of 
our LM data, reinforcing the conclusion that this system 
does not present any freezing transition at A = oo. 

The pair correlations exhibit the same sort of behav¬ 
ior as in the 2d case: only the 1st few nearest neighbors 




A 


FIG. 7: Spin glass susceptibility (top) and correlation length 
(bottom) as computed from the LM approach (lines) or mea¬ 
sured in the MC simulations (points), for the LGF interaction 
on the cubic lattice. 


tend to be strongly antiferromagnetically correlated, but 
no correlations develop at large distances as the coupling 
strength is increased, and the system remains paramag¬ 
netic. 


V. SPECTRAL PROPERTIES 

The A“^ = 0 transition can be considered from the 
point of view of the interaction matrix and CZD, 

as an example of euclidean random matrix (ERM) un¬ 
like the traditional random matrices, where different en¬ 
tries of the matrix are uncorrelated, ERM’s are defined 
by a function of the distance between two points /(c), 
where the randomness in the entries is induced by the 
randomness of the underlying point pattern {ri}. These 
random matrices have been studied for certain classes of 
functions /— , and some classical results are available. 
Our degree of understanding of this subject is not com¬ 
parable to that of the classical (e.g. GOE,GUE, Wishart) 
ensembles21 with most results coming from exact diago- 
nalisation and approximation s . 
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FIG. 8: Scaling of the number of vanishing eigenvalues of the 
matrix B defined on the text and of the spin glass susceptibil¬ 
ity (inset) with the number of particles for the LGF in three 
dimensions at 1/A = 0. 



X 

FIG. 9: Ground state eigenvector showing a trimer for a par¬ 
ticular disorder realisation using the LGF as interaction on a 
lattice of size L — 32 with N = 102 particles. The compo¬ 
nents of the eigenvector are proportional to the radius of the 
circles, which are centered on the corresponding spin position. 
Red (black) sign indicates a positive (negative) sign. 


Unfortunately due to the long-range nature of the log- 
interaction, many of the methods to analyse the spectral 
properties presented in Ref. [T^ do not apply directly to 
our case. However, a phenomenological picture of the 
low- and high-lying eigenstates of the matrix Jij can be 
established transparently. 

Let us start from the large positive eigenvalues. Since 
Jij is constant in sign, the Frobenius-Perron theorem 
states that a highest eigenvector is nodeless. To a rea¬ 
sonable approximation, it is fully delocalised, 

~ (l/y]v,...,l/ViV). (31) 

The associated eigenvalue is 

N 

Amax^ylniV. (32) 


with an inverse participation ratio of 1/A^. 

The second-to-highest eigenvalue is also associated to 
a delocalised eigenvector, which is now a wave with wave¬ 
length < L. At these length scales the randomness of the 
point process plays little role. A finite fraction (possi¬ 
bly all) of the eigenstates containing the largest eigenval¬ 
ues are delocalised, they correspond to long-wavelength 
charge-density variations. The average spectral density 
(DOS) of the LGF (fT6l) interaction matrices is shown on 
top panels of Figs. ITUlandfTTl in the limits of high {x = 
0.125) and low density (‘continuum limit’, x = 2“^^), 
respectively. 

Guided by the numerics, we see that the eigenvectors 
corresponding to the most negative eigenvalues are lo¬ 
calised eigenvectors: most of the weight is concentrated 
on 0(1) spins. This leads us to consider isolated perco¬ 
lation animals. 

The simplest (and, for small x, the most abundant) of 
these is the dimer. A well-isolated dimer supports two 
eigenvalues: an antisymmetric and a symmetric one. The 
antisymmetric one, 

<^(0) = (i/^,_l/v/2,0,...,0) (33) 

corresponds to the smallest eigenvalue. In fact, since the 
closest pair is located one lattice spacing away J 12 ^ In L 
and the lowest eigenvalue is 

A„,i„ ~ - ln(L) + 0(1) ~ i \n{N/L^) - ^ ln(N) + 0(1). 

(34) 

At fixed density, N/L?, the lowest eigenvalue depends 
logarithmically on the system size. 

For a well isolated dimer, say at distance r from the 
closest spin, the effect of neglecting the rest of the spins 
appears as a correction 0(l/r). 

We now consider how big this isolation distance r is. 
By the usual arguments of percolation theory, one can 
estimate the expected number of isolated dimers as 

nd{r) = L^2x'^{l-xY''^, (35) 


where we have approximated the number of lattice sites 
in a circle of size r with 7rr^. Therefore the most iso¬ 
lated dimer (the solution of the equation nd{r) = 1) is 
surrounded by an empty area of size 


r[L) 


^2 \n{xLY2) 

\/7rln(l/(l-a;)) 


(36) 


Note the extremely slow dependence r[L) ^ Vln L. 

Inserting L = 120 and x = 0.1, which is about the 
largest sizes considered in our numerics, r = 4.13, which 
can hardly be called isolated. 

The isolation effect would be much more pronounced 
for X = = 1,200, for which r = 18.4. Other¬ 

wise, one needs to consider the ground states of more 
complicated lattice animals, like trimers, snakes, squares 
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I + 0.5ln(N) 

FIG. 10: Spectral density (top) and average Y (Eq. 1371 bot¬ 
tom) for a fraction of a; = 0.125 occupied sites in the lattice, 
using Jij as defined in mi), the LGF interaction. The inset 
shows the fluctuations of Y. 



FIG. 11: Spectral density (top) and average Y (bottom) for 
a fraction of a; = 2“^^ occupied sites in the lattice, using Jij 
as defined in mi, the LGF interaction. The inset shows the 
fluctuations of Y. 


etc. As an example, a ground state eigenvectors for one 
disorder realisation is shown on Fig. [5] 

This problem becomes quickly analytically prohibitive. 
However the fact that the ground state is localised on 
some lattice animal appears robust: on the graphs we 
consider, the smallest eigenvalue is ^ — InL and the IPR 
is 0(1). 

With the lower end of the spectrum localised and the 
high-end delocalised, it is a natural question whether 
there exists a mobility edge separating the two limits. 
In order to study the transition we have looked at the 
inverse participation ratio as a function of the eigenvalue 
A: 


Y{X) 


IPR„=^ui 

i 


1 

pW 


^IPR„5(A-A„), 


(37) 


where Aa and Vai are eigenvalues and normalized eigen¬ 
vectors of Jij respectively. We consider the average 
[T](A) and fluctuations a{Y){X).— A mobility edge would 
be signaled by the divergence of the fluctuations of Y (A) 
at a certain Ac. Numerical diagonalization of Jij does not 
indicate such a transition: the two limits appear to be 


separated by a crossover. The bottom panels on Figs. fTOl 
and m show, respectively for a high and low density of 
particles, the average Y{X), while the insets display the 
fluctuations of Y (A). The spectral properties of the LGF 
in d = 3 turn out to be very similar to the d = 2 case 
(not shown). 

A detailed study of this ERM ensemble would be de¬ 
sirable and is left for future work. 


VI. PAIR CORRELATIONS AND SCREENING 
A. Analytical theory of screening 


Away from the T —^ 0 limit of the microscopic model, 
excitations of the non-orphan tetrahedra out of their mo¬ 
mentless state carry gauge charge, which leads to a vari¬ 
ant of Debye screening, with the special feature that the 
gaplessness of the charge excitations leads to a some¬ 
what unusual temperature dependence of the screening 
length^. 

In addition to this, even in the limit T —>■ 0 studied 
here, we encounter an additional type of screening. This 
occurs on account of the long-range uniformly antiferro- 
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1 / A 

FIG. 12: Uniform susceptibility as computed from MC sim¬ 
ulations (dots) compared to a weak coupling expansion av¬ 
eraged over 200 disorder realisations (dashed lines). The 
Curie-Weiss constant increases approximately linearly with 
the number of particles (inset), yielding a vanishing radius of 
convergence of the weak coupling expansion already at leading 
order. 


magnetic Coulomb interaction between the orphan spins, 
whose existence is the distinguishing property of the ran¬ 
dom Coulomb antiferromagnet. It again exhibits a Debye 
form, although distinct from the setting of mobile charges 
in which Debye screening is normally considered, as here 
it is the (continuous) flavour of the charges - the orien¬ 
tation of the orphan spin whose orientation is free but 
whose location is fixed - which is the dynamical degree 
of freedom. 

This can be seen directly in a weak-coupling expan¬ 
sion, which in Coulomb systems has a vanishing radius 
of convergence in the thermodynamic limit, as is easily 
verified in our simulations. Fig. 1121 

To elucidate the role of screening, we compute the dis¬ 
order averaged correlator between two spins at Tq and . 
Consider the Hamiltonian 

77 = — ^ ^ Jijni ■ rij, (38) 


(with a factor of 3 to guarantee that = 

3/3= 1). Thus 


-aJ 

Cab = Sab-{a\-j-^\b) (41) 

where we use a matrix notation (a| J|&) = Jab- For sim¬ 
plicity we will not write the Sab term, which only affects 
the result for the self-correlation (it will return to be 
important when we discuss the LM approximation again 
later). The correlation function between a and b depends 
also on the positions of all the other points X 2 , ■■■,xn so 
it should be written as C{xa,xb\x 2 , ---^xn)- 

This Gaussian approximation is equivalent to the re¬ 
summation of a set of diagrams in which there are no 
internal loops, dubbed “chain diagrams.” This approxi¬ 
mation is justified in the limit of small a, in which spins 
are rarely polarized along some direction and the hard- 
spin constraint is not so important. 

This result holds for each disorder realization. We now 
take the average over realizations (leaving the question 
of whether this is representative of the distribution or 
not for later) keeping fixed the position of the two spins 
a, b. For doing this, it is convenient to go back to the 
geometric expansions and define 



d^-'^x 

gN-2 C{Xa,Xb\xi,...,XN-2) 


(42) 


where Xi are the locations of the other TV — 2 spins and 
S = LJ . We have relaxed the constraint that points be 
located on a square lattice, which is immaterial in our 
high temperature, low-dilution expansion. 

Unfortunately it is difficult to see what the distribution 
of J induced by the random positions is, but we can 
expand the Gaussian result in powers of a and do the 
average term by term. 

We get 


E [Cab] 


i 
3 

^ ^ E \JaiJij Jjb\ 



j E 


(43) 


where Jij are given by either the Log or the LGF and 
we will eventually set a = 1. The correlation function 
between two spins, for fixed disorder is: 


Cab — ij^a ' ‘bib) — 
= i <5(1 - nf)(na • 


As it is not the hard spin constraint which is central 
to the physics of screening, we substitute it with some¬ 
thing more manageable (analogously to the LM method, 
but without imposing self-consistency). Representing the 
delta function with a Gaussian term 


(27r/3)3/2® 


Now, term by term we obtain objects like 

f dJ'x 

E ^JaiJib = {N-2) / —J{Xa-x)J{x-Xb) 

. i J 

— pj d^xJ{Xa — x)J{x — Xb) (44) 

where p = (TV — 2)/S ~ N/S is the density of points. 
Fourier transforming, 

p f d^xJ{Xa — x)J{x — Xb) 

— n f j 2 d q d q T iq{x^-x)+iq'{x-Xb)( AC-\ 
“ ^7 (27r)2(27r)2'^«'^« ® 

[ ^ d j2 iq{xa.-Xb) 

J (27r)2''^" 




(40) 


P 


(46) 
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The geometric series obtained thus for E [Cab] yields 


E [Cab] 


^ibCa-xh) (a/3)Jg 
(27r)2 1 + {ap/3)Jq 


(47) 


Now, for both Log and the LGF, Jg ~ c/q^ (c is a con¬ 
stant of 0(1))^ so that at small a we have approximately 


E [Cab] 


(Pq ig(x„-x^) (ca/3) 

{2TTy q^ + {cap/3)' 


(48) 


This leads to 

E [Cab] ^ {-2a/3)Ko{r^/cap/3) (49) 

which exhibits a screening length 

e = l/v'cap/3. (50) 

As both a and c are 0(1) this shows (not surprisingly) 
that the screening length is proportional to the 1/^^- 



r 


FIG. 13: Correlation function exhibiting screening: numerical 
results (for a single disorder realization with N = 200 points 
on a square of unit size) compared to the predicted analytical 
form from the chain diagrams. 


Note that in this approximation, for & <C C the cor¬ 
relation function C{r) ^ 1, which is not physical for unit 
length spins. This is an artefact resulting from substitut¬ 
ing the hard spin constraint with a quadratic confining 
potential. Therefore this approximation is internally con¬ 
sistent only for Ca^b ^ Cj where it predicts an exponential 
damping of the correlations but we note that the large 
anticorrelations at short distance due to strongly coupled 
spins close to one another put these into a state with van¬ 
ishing total spin, which - physically correctly - screens 
their joint field at larger distances. 


B. A random scattering picture 


To gain some insight into this, we develop an analogy 
with wave propagation in disordered media, which sug¬ 
gests that no transition exists. The basic observation is 
that the interaction is simply related to the inverse of the 
Laplacian, the propagator of a free particle on the lattice: 

Considering that 

Jij = (51) 

properly regularized (particularly important is the con¬ 
dition that Jii = 0), we can rewrite the expression (HTl) 
as 

Cab = Sab-^{a] _^^^^_^ ]b), (52) 

where E = 0 and 

i 

is a random potential. This can be established by ex¬ 
panding in powers of a. 

Thus C is (proportional to) the propagator for a wave 
in a two-dimensional box with randomly placed point-like 
scatterers^2iHj at energy E = 0. 

The precise form of the mapping is the following: the 
correlation function 

3 

- {naUb), (54) 

a 

is the amplitude of a signal sent from the scatterer a to 
the scatterer 6, considering all order processes bouncing 
over all the N scatterers. In case a = b the direct path 
from a to 6 needs to be neglected. This is a form of 
renormalization of the scattering problem which is always 
necessary in the point-like (or s-wave) scattering limit^i. 

Once the renormalization procedure is done, the prob¬ 
lem we are left with corresponds to the propagation of 
a scalar wave, damped by a scattering section for every 
typical realization of disorder. Without repeating the 
classical treatment of this phenomenon we can say that 
the signals (spin-spin correlations) must be screened for 
any a, the screening length (measured in units of 1/^/p) 
being a decreasing function of a. Even if not precisely of 
the form (ISUll for small-a, it seems to diverge like 1/^/a. 

This is valid both for the coherent field E [Cob] and the 
incoherent field E [C^j,] — E [Cob] , although the scatter¬ 
ing sections (and hence the damping/correlation lengths) 
might have different values. This analogy makes us re¬ 
alize that in this approximation there is no transition 
irrespective of the value of a or p, and this is consistent 
with numerical results. 

This analogy extends also to the LM limit. Consid¬ 
ering a small-a series expansion for the spin correlation 
function: 


The final question we address concerns the fluctuations 
of the random quantity (HTl) and whether these may sig¬ 
nal any phase transition even when the mean does not. 


haCabhb — <^ab^a aJabEa Jai , Jib a Jai , Jij , Jjb~\~--- 

I ^2 ’ ^2 ! L q 

( 55 ) 
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(recall that in LM a is scaled by a factor 1/m, hence 
the factor of 3 of the previous paragraphs is absent here) 
where the extra factors of ha need to be chosen in such 
a way that 

Qi = (riini) = 1. (56) 

Cab is then proportional to the propagator 

Gab = ^ab-a(a|—^|6), (57) 

where h is the diagonal matrix with diagonal entries 
E = 0 and 

V{x) = a'^^5{x - Xi), (58) 

i 

where the renormalized value = 0 is intended. 

This is a scattering problem over point-like scatterers, 
where now each scatterer has different scattering ampli¬ 
tude. This modification should not change the physi¬ 
cal analogy of the problem. This is again a scatter¬ 
ing problem of a scalar wave over point-like scatterers. 
The propagation of the wave is attenuated over distance 
in the usual exponential fashion. Therefore, if a phase 
transition exists, it is not mirrored in the divergence of 
the correlation length. Conversely, as this treatment is 
closely related to the LM one (rather than the Heisenberg 
model), on account of the softening of the hard constraint 
to a Gaussian one, we would not expect a transition at 
finite value of a. 


VII. DISCUSSION 

We have studied the effective theory describing disor¬ 
der in the form of quenched non-magnetic impurities, in 
the topological Coulomb phase, on a lattice with bipar¬ 
tite dual. Interactions in the effective picture are long- 
ranged, and to the best of our knowledge this is the first 
study available of such a model. 

A. A freezing transition? 

Our results show that any freezing transition, if it ex¬ 
ists, is extremely tenuous. In d = 2, for LM there does 
not appear to be freezing for any finite coupling, with 
a nice scaling collapse of the data at A = oo indicating 
freezing to take place in this limit. 

The relation of this result to a finite number of spin 
components is the following. Firstly, our Heisenberg sim¬ 
ulations cannot access a freezing transition, but they do 
show a greater tendency towards glassiness than LM, 
with both a larger spin-glass correlation length and an 
enhanced tendency for the curves to cross. 

This is in keeping with the general expectation^! for 
the more constrained Heisenberg model to freeze before 


the soft spins do (anf after an Ising model might). If there 
is a freezing transition at < oo, it will still be at phan- 
tastically large coupling Ac > 100. The delicate nature 
of all of these phenomena is further underscored by the 
dependence on finite-size choices, which may lead to an 
entirely different set of instabilities. Similarly, the analyt¬ 
ical approaches, in particular the mapping to a quantum 
scattering problem, see little indication of a transition. 

The tendency towards freezing seems to be even weaker 
in d = 3, perhaps surprisingly so, given the freezing tran¬ 
sition is more robust in higher dimension for the instances 
of canonical spin glasses. However, unlike in these cases, 
our distribution of the intersite couplings is dimensional¬ 
ity dependent, and in particular becomes ’shorter-ranged’ 
as the power law of the decay of the Coulomb law grows 
with d (while, of course, the power law with which the 
number of distant spins grows, increases). 

The weak tendency towards freezing is in keeping with 
the fact that our model is not easily deformed into one 
of the standard spin glass models. On one hand, in¬ 
creasing the range of the interaction towards the extreme 
of doing away with any notion of distance and assign¬ 
ing equal coupling between all the spins yields simply a 
global charge-neutrality constraint (which, at any rate, is 
already enforced microscopically) and therefore preserves 
a microcanonical version of a perfect paramagnet. If the 
coupling is restricted to nearest-neighbours only, we in¬ 
stead get a combination of percolation physics and that of 
the standard Neel state for a bipartite antiferromagnets, 
where any tendency towards disorder is a dimensionality 
effect, and glassiness is nowhere to be seen. 

The tendency towards glassiness is therefore necessar¬ 
ily due to a combination of the non-constancy of the log¬ 
arithmic interaction - which, helpfully, is not bounded 
as r oo, along with its long range. Studying mod¬ 
els exhibiting this pair of ingredients more systemati¬ 
cally is surely an interesting avenue for future research. 
We would like to emphasize, in particular, that the phe¬ 
nomenon of screening we have discussed has no counter¬ 
part in the literature on conventional spin glasses, where 
the random choice of the sign of the interactions does not 
allow the identification of an underlying charge structure. 

In this sense, our model is much closer to those fa¬ 
miliar from the study of Coulomb glasses, although the 
differences here are again considerable. We have vector 
charges rather than Ising (positive or negative) ones; dis¬ 
order appears in the form of random but fixed locations 
rather than fixed on-site potentials for charges not bound 
to a particular site. It is intriguing that such a variation 
of a classic Coulomb glass appears entirely naturally in 
frustrated magnetism. 


B. Freezing in frustrated magnetic materials 

With Heisenberg spins placed at random sites of the 
pyrochlore-slab lattice (also known as the SCGO lattice) 
and a particular, microscopically determined value of A, 
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the d = 2 case of our Coulomb antiferromagnet corre¬ 
sponds, up to the sublattice-dependent inversion factor 
mentioned earlier, to the T —>■ 0 limit of the physics of 
orphan-spins created when a pair of Ga impurities substi¬ 
tutes for two of the three Cr spins in a triangular simplex 
of this lattice. Although experimental interest in SCGO 
dates back to the 80s and played a key role in stimulat¬ 
ing experimental and theoretical interest in the area of 
highly frustrated magnetismii, the behaviour of SCGO 
is reasonably well-understood in theoretical terms only in 
the broad Coulomb spin-liquid regime down to about a 
hundredth of the exchange energy scale (of order 500K) 
between the Cr spins. The magnetic response in this 
regime can be modeled in a rather detailed way as being 
made up as the response of a pure Coulomb spin-liquid 
superposed with the Curie-tails associated with vacancy- 
induced “orphan-spin” degrees of freedom^^— that carry 
an effective fractional spin^i^ and leave their imprint 
on NMR lineshapesi^ and bulk-susceptibility^>^‘i^ in the 
Coulomb spin-liquid phase. In contrast, the physics at 
very low temperatures (of order 5K or lower) is still not 
very well understood, with intriguing but largely unex¬ 
plained reports of observed glassy behaviour even at very 
low densities of Ga impuritiesi^d^, which appears to in¬ 
volve only the freezing of a fraction of its degrees of free¬ 
dom. 

Our model retains the key feature of the T —>■ 0 limit of 
the effective model, namely the long-range Coulomb form 
of the effective exchange couplings, but does not retain 
the detailed geometry of these orphan-spins in SCGO, ex¬ 
cept for the sublattice-dependent inversion that connects 
the degrees of freedom of our Coulomb antiferromagnet 
with the underlying physics of these orphan-spins. 

Bearing all this in mind, the usual caveat about ide¬ 
alised models for frustrated systems applies to our study 
as well: Our starting Hamiltonian of a classical nearest- 
neighbour Heisenberg model does not include a num¬ 
ber of aspects - further-neighbour interactions, single¬ 
ion anisotropies, non-commutation of spin components - 
all of which give rise to interesting, generally non-glassy, 
physics of their own. If and when these energy scales 
dominate over our the instabilities of the idealised model, 
it is the former which will likely show up more promi¬ 
nently in experiment. 

In addition, in our case, the critical coupling Ac, even 
if it is not infinite, is hard to attain in any microscopic 
model. Indeed, for the checkerboard lattice, one obtains 
A = I/47r from a microscopic calculation, easily within a 
very short-range correlated regime. 

At any finite temperature, which is all that can be ac¬ 
cessed experimentally for the time being, the Coulomb 
interactions obtain a finite-screening length due to the 
thermal excitation of charges even in non-orphan tetrahe- 
dra. Following the general lore on spin freezing, this pre¬ 
cludes even canonical Heisenberg spin glassiness. For this 
reason, the abovementioned A-independence of a freez¬ 
ing transition in d = 2 is not going to carry over directly 
to the experimental compound. 




FIG. 14: (Left)The sign function on the 1st quadrant for a 
lattice of side L = 100 used to resum the correlations. (Right) 
The scaling of the new susceptibility proposed to describe the 
ordering occuring with the Log interaction. 



FIG. 15: Numerically obtained Fourier transform of Log (red) 
and LGF (blue) for a lattice of size L = 100. Circles connected 
by lines indicate the edge ky = 0, while squares indicate the 
diagonal kx = ky. The index n labeling the x axis indicates 
the index of the wave vector: kx = ^n. The inset shows in 
more detail the region near the global minimum. 


However, real systems will only be quasi-2(i, with resid¬ 
ual couplings between the two-dimensional layers. In¬ 
deed, for the case of SCGO, dilution also breaks up the 
tightly bound singlets of the dimers of Cr ions which iso¬ 
late the kagome-triangle-kagome trilayers from one an¬ 
other. The consequences of coupling in the third dimen¬ 
sion remain an interesting yet completely open topic for 
future study. 


C. Connection to other models 

More broadly, perhaps the most pleasing aspect of this 
work is how it naturally connects (with) a number of de¬ 
formations of well-known problems-the scattering prob¬ 
lem, Coulomb glass physics, or random matrix theory. 
In particular, we have identified a straightforward way 
of obtaining a Euclidean random matrix problem from 
a simple magnetic model where long-range interactions 
emerge naturally. We hope that this will motivate fur¬ 
ther work on any (and perhaps all) of these problems. 
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Appendix B: Fully Occupied Lattice 


Proposition 1 in Ref. [2^ states that if J{k) is the 
Fourier transform of the interaction matrix J, then a 
minimizer /cq for J{k) determines a modulated ground 
state for the system with that wavevector. 

The Fourier transform of the LGF at nonzero wavevec¬ 
tor is readily read from its definition, Eq. (HU): 


jhGF{k) 


1 

2 — COs(fca;) — COs{ky) 


(Bl) 


FIG. 16: The average on each binning block of the spin glass 
susceptibilities plotted against the logarithm (base 2) of the 
size of the corresponding binning block. Data shown here 
corresponds to MC simulation of the LGF in the cubic lattice 
at A = 200. 
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Appendix A: Non-Glassiness for the Log Interaction 


The pair correlation profiles for the Log interaction 
exhibit a structure hinting on the way pair correlations 
should be summed in order to define a generalized suscep¬ 
tibility describing the order present on this system. This 
order reflects the symmetry of the interaction, which is 
anisotropic, but has the symmetries of the square lattice. 

We define a sign function, 9{x, y), which on each quad¬ 
rant has alternating values ±1 on suscessive “square 
frames” of fixed width of 2 lattice sites for any L. As¬ 
suming (x, y) on the 1st quadrant, this function has the 
profile pictured on the left panel of Fig. [TH 
The corresponding susceptibility reads: 


X = 




(Al) 


Square brackets denote as usual disorder average. This 
susceptibility diverges with system size, and its scaling in 
MG simulations is shown on the right panel of Fig. [TT] 
the same behavior is found in the LM. 


which has a minimum at fc = (tt, tt), thence we find “con¬ 
ventional” antiferromagnetic order. 

For the Log interaction, we are not able to find an 
analytical expression for its Fourier transform, but nu¬ 
merical results show that the global minima happen at 
k = (tt, 0) or (0, tt), which explains the striped phase for 
the fully occupied lattice. The non-analyticity of the dis¬ 
tance function periodized by the functions min(x, L — x) 
or min(i/, L — y) (which is seen as a discontinuity in the 
derivative along the lines x = L oy y = L j^') gives rise 
to “ringing” in J^ogik), a line of alternating local max¬ 
ima and minima appear along k^ = 0 or ky = 0. The 
new global minimum is shifted from (tt, tt) to the edges 
of these lines, as shown in Fig. [TS] 


Appendix C: Verifying Equilibration 


Our simulations require exploring a region of very 
high coupling, A. In this case it is important to en¬ 
sure that equilibrium is attained. To test this, we bin 
the data for the spin glass susceptibility. This bin¬ 
ning consists of subdividing the total number of mea¬ 
surements, Nm, in contiguous bins of successive sizes: 
1,1, 2,4,8,..., iVm/4, Vm/2. The average for each bin 
is then plotted against the logarithm of the bin size 
(Fig. [1^1). Equilibrium is diagnosed by at least the last 
3 bin averages agreeing within the interval set by their 
error bars. 

The final equilibrium values used consist of the average 
of the last half of the measurements made in the simula¬ 
tion, Nm/‘2. 
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